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EFFECTS  OF  ANNEALING  ON  THE 
ISOCONFIGUBATIONAL  FLOW  OF  A  METALLIC  GLASS 


S.  S.  Tsao*  and  F.  Spaepen 


Division  of  Applied  Sciences,  Harvard  Univeraity 
29  Oxford  Street.  Caabridge,  MA  02138,  C.S.A. 


ABSTRACT 


The  i soconf  ignrati onal  viscosity  of  amorphous  Pd^  ^cijSi^gi  in 
various  states  of  strnctaral  relaxation  was  investigated  over  a  wide 
temperature  range.  It  was  found  to  be  thermal  history  dependent  and 
non-Arrhenian.  The  lower  the  glass  fictive  temperature,  and  hence 
the  more  ordered  the  glass,  the  greater  the  departure  from  Arrheuian 
behavior.  The  non-Arrhenian  behavior  shows  that  the  temperature 
dependence  of  the  isoconf igurational  viscosity  cannot  be  ascribed 
entirely  to  the  temperature  dependence  of  the  jump  frequency  of  the 
flow  defect.  The  observed  crossover  of  isoconf igurational  curves 
obtained  from  different  structural  states  implies  that  the 
isoconf igurational  viscosity  cannot  be  a  function  of  only  a  single 
structural  parameter,  such  as  the  average  free  volume  of  the 
specimen. 


•Present  address:  Sandia  National  Laboratories,  Albuquerque,  NM  87115,  U.S.A. 


1. 


INTRODUCTION 


For  measurements  of  the  temperature  dependence  of  the  viscosity,  q(T).  of 
an  amorphous  metal  in  its  (unstable)  glassy  state  to  be  physically  meaningful, 
they  must  be  isoconf i aura t ions 1 . 

For  example,  if  the  viscosity  were  simply  measured  at  increasingly  higher 
temperatures,  as  a  result  of  structural  relaxation,  which  causes  continuous 
viscosity  increase,  the  apparent  activation  energy  could  be  too  low,  and 
physically  meaningless. 

In  an  ideal  isoconf igurational  flow  experiment,  the  viscosity  of  the  same 
structure  (i.e.,  one  in  which  the  relative  positions  of  the  atoms  remain,  at 
least  statistically,  the  same)  is  measured  as  a  function  of  temperature.  In 
the  Newtonian  regime,  the  viscosity  can  be  written  as  [1]: 


(n^:  concentration  of  atomic  scale  sites  producing  shear,  or  ''flow  defects''; 
bf':  stress  unbiased  jump  frequency  of  a  shear  defect;  r0:  shear  strain  per 
jump  of  a  defect;  vo;  volume  of  a  defect).  The  temperature  dependence  of  q 
is,  in  general,  dominated  by  the  product  n^k^'.  The  atomistic  interpretation 
of  the  results  from  an  ideal  isoconf igurational  experiment  is  then 
straightforward:  the  number  of  defects  remains  constant;  therefore  the 
temperature  dependence  of  the  isoconf igurational  viscosity  nitQ(T)  reflects 
only  changes  in  the  jump  frequency  k^’. 

In  practice,  however,  it  is  difficult  to  ascertain  experimentally  that 
these  ideal  isoconf igurational  conditions  are  achieved.  All 
' 'isoconf igurational ' '  measurements  are  therefore  only  defined  in  an 


operational  sense;  i.e.,  the  effects  of  irreversible  strnctnral  relaxation  are 
eliminated,  bnt  nothing  is  a  priori  known  about  possible  reversible  changes  in 
the  atomic  structure  of  the  glass. 

In  one  type  of  isoconf igurational  experiment,  the  temperature  of  the 
system,  initially  held  at  equilibrium,  is  rapidly  changed  by  a  small  amount, 
and  the  instantaneous  viscosity  change  determined  by  extrapolation  of  the 
viscosity  transient  back  to  the  time  when  the  new  temperature  was  reached  [2] . 
This  type  of  isoconf igurational  experiment  is  suitable  for  easy  glass  formers, 
such  as  silicates  and  polymers. 

For  metallic  glasses,  this  technique  is  less  suitable,  since  the  approach 
to  (metastable)  equilibrium  is  often  accompanied  by  rapid  crystallization. 
Isoconf igurational  experiments  on  metallic  glasses  must  therefore  be  performed 
mostly  below  the  glass  transition  temperature,  T^.  Although  structural 
relaxation  in  these  experiments  leads  to  viscosity  increases  with  time  (at 
temperatures  below  the  fictive  temperature,  T^) ,  and  even  linearly  far  from 
equilibrium,  the  relative  rate  of  increase  d(ln  q)/dt,  decreases  with  time. 

By  structurally  relaxing  the  specimen  for  a  sufficiently  long  time,  one  can 
measure  the  viscosity  at  a  number  of  different  temperatures,  return  to  the 

original  temperature  and  find  that  the  viscosity  has  not  changed  measurably. 

-4 

Since  a  typical  viscosity  measurement  requires  a  total  strain  of  only  10  , 
the  total  number  of  atomic  jumps  during  the  experiment  (assuming  y  ~  1)  is 
only  one  in  10*.  This  explains  why  isoconf igurational  measurements  of  this 
type  are  possible,  since  the  effects  of  irreversible  atomic  jumps  on  the 
structure  can  be  minimized.  The  basic  feature  of  this  type  of 
isoconf igurational  experiment  is  thus  reproducibility  of  the  measurements  on 
temperature  cycling. 

Several  investigators  have  performed  the  second  type  of  experiment  on 


metallic  glass  specimens  far  from  equilibrium  and  found  an  Arrhenian 
temperature  dependence  of  Taub  and  Spaepen  measured  activation 

energies,  Q^s0>  of  192kJ/mole  and  240kJ/mole  in  PdgjSijg  (31  and 
Pd^^  jCugSijg  j  [41,  respectively.  Mulder  et  al.  measured  Q^#0  of  250kJ/mole 

in  Ni36Fe32Cr14P12B6  ^  and  Fe40Ni40B20 

As  discussed  in  an  earlier  paper  [71,  the  straightforward  interpretation 

of  Q^s0  as  simply  reflecting  Arrhenian  temperature  dependence  of  the  jump 

frequency  k^'  (i.e.  n^  constant)  poses  a  problem.  In  that  case,  the 

isoconf igurational  viscosity  would  also  have  an  Arrhenian  form: 

\so(T)  -  %°  exp  (Qi,o/kT)  (2> 

with  q*  =  kT/nf\)D(yovo)  only  weakly  temperature  dependent;  ;  Debye 
frequency. 

For  the  first  (i.e.,  least  relaxed)  set  of  the  Taub  and  Spaepen 

isoconf  igurational  data  on  Pdg2Sijg  [31,  q*S0  “  1 .5xlO-^Nsm-2 ,  which  is  about 

—3  —2 

four  orders  of  magnitude  lower  than  the  expected  value  of  about  10  Nsm 
(viscosity  of  a  liquid  at  high  temperatures),  and  furthermore,  is  more  than 
two  orders  of  magnitude  lower  than  the  absolute  minimum  value  that  can  be 
obtained  from  equation  (2)  (using  n^  =  1/0,  \)g  =  10*2sec  * ,  Y0v0  *  3.40  [81; 

Q:  atomic  volume) . 

In  view  of  this  large  discrepancy  in  the  value  of  the  pre-exponential,  it 
seemed  that  iaoconf igurational  flow  could  not  be  considered  to  be  a  singly 
activated  Arrhenian  process.  To  teat  this  we  performed  isoconf igurational 
experimenta  over  a  wider  temperature  range  and  closer  to  T^.  In  the 
experiments  of  Taub  and  Spaepen  [3,4]  and  Mulder  et  al.  15,6],  each  apparently 
Arrhenian  isoconf igurational  data  set  only  spanned  401  or  less. 
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These  investigators  also  found  to  be  tbe  sane  for  a  noaber  of 

stroctnral  states,  i.e.  independent  of  thermal  history.  In  this  paper,  ve 
report  results  of  experiments  designed  to  investigate  the  iaoconf ignrational 
behavior  over  a  wide  temperature  range,  particularly  at  temperatures  close  to 
Tg,  end  to  determine  whether  Q^so  remains  independent  of  thermal  history  if 
the  system  is  subjected  to  more  markedly  different  annealing  treatments. 

2.  EXPERIMENTAL  PROCEDURE 

Viscosity  measurements  were  made  using  a  tensile  creep  apparatus  [9], 
modified  as  described  previously  [10].  The  specimens  were  P d^  - Cu^ S ij ^  5 
disk-quenched  ribbons  and  wires.  The  specimens  were  under  load  during  the 
preanneals  necessary  to  stabilize  the  amorphous  structure  prior  to 
isoconf igurational  testing.  Stresses  were  kept  low,  from  21-144MPa  to  ensure 
Newtonian  viscous  flow.  Due  to  the  strong  temperature  dependence  of  the 
strain  rate,  test  times  varied  from  a  few  minutes  at  high  temperatures  to  a 
few  days  at  low  temperatures.  In  early  experiments,  the  stress  was  increased 
or  decreased  during  the  temperature  cycling,  depending  on  the  direction  of  the 
temperature  change,  in  order  to  extend  the  viscosity  range  of  measurements 
within  one  isoconf igurational  set. 

In  the  last  sets  of  experiments,  »iigo(T)  was  better  characterized  and 
wide-ranging  i soconf igurational  data  sets  (over  90K,  and  three  orders  of 
magnitude  in  (T))  were  measured  without  any  load  changes.  It  is  desirable 
to  maintain  the  same  load  during  temperature  cycling  since  a  load  change 
induces  anelastic  strains  that  must  be  taken  into  account.  For  example,  from 
loading-unloading  experiments,  we  found  that  the  anelastic  strain  contributes 
measurably  to  flow  for  about  a  day  at  527K,  and  for  about  6  minutea  at  615K. 


Thus  longer  test  times  ere  required  for  measurements  after  a  load  change,  in 
order  to  allow  for  anelastic  transient  decay.  Shorter  test  times  are 
desirable  to  avoid  structural  relaxation. 

3 .  RESULTS 

Figure  1  summarizes  the  results  on  a  Pd^^  ^Cu^Si^  j  ribbon  first 
annealed  at  529K  for  457  hours.  Point  1  shows  the  final  viscosity.  The 
numbers  indicate  the  order  of  the  tests.  After  the  measurement  at  point  4 
with  a  test  time  of  24  hoars,  the  shear  stress  was  decreased  from  144MPa  to 
96MPa  and  the  sample  held  at  the  same  low  temperature  of  501K  for  an 
additional  97  hours  to  allow  for  anelastic  transient  decay.  Point  5  and,  to  a 
lesser  degree  point  6,  are  too  high  due  to  remaining  transients  from  the  load 
decrement  at  a  temperature  (501K)  lower  than  the  stabilization  temperature 
(test  1,  529K) . 

Complete  transient  decay  has  likely  occurred  prior  to  the  measurement  at 
point  7.  Points  1-20  trace  out  a  first  isoconf igurational  curve,  as 
demonstrated  by  cycling  the  temperature  and  reproducing  the  viscosity  values. 
After  the  viscosity  measurement  at  point  20,  the  specimen  was  allowed  to  relax 
at  588K  for  40  hours.  Point  21  shows  the  final  viscosity.  A  second 
isoconf igurational  curve  (points  21-30)  was  then  determined. 

The  results  on  a  wire  specimen,  preannealed  under  load  for  73  hours  at 
593K,  are  summarized  in  Fig.  2.  Point  4  is  too  low  due  to  transients  after 
increasing  the  shear  stress  from  22MPa  to  67MPa. 

Figure  3  shows  the  results  on  a  ribbon  specimen  first  annealed  at  593K 
for  4  hours.  Point  1  is  the  viscosity  at  the  end  of  the  preanneal.  Points 
1-4  determine  a  first  isoconf igurational  set  A,  points  7-17  belong  to  a  second 


set  B.  The  viscosity  value  shown  at  point  17  was  measnred  at  the  start  of  a 

19-hour  anneal  at  593K.  Point  18  shows  the  final  value.  Points  18-30  trace 

out  a  third  isoconf igurational  curve  C.  The  value  shown  at  point  30  was 

measured  at  the  start  of  a  24-minute  anneal  at  615K.  Since  this  value  is 

greater  than  the  equilibrium  viscosity,  q  ,  at  this  temperature  [10],  the 

eq 

viscosity  decreased  with  time  to  that  shown  at  point  31.  A  fourth 

isoconf igurational  set  D  (points  31-40)  was  then  determined.  The  same 

specimen  was  then  once  again  annealed  at  593K  for  24  hours.  Points  42  and  43 

show  the  viscosity  values  at  the  start  and  end  of  the  anneal,  respectively.  A 

fifth  isoconf igurational  set  E  (points  43-52)  was  determined  before  the 

specimen  was  cooled  to  room  temperature.  The  shear  stress  on  the  specimen  was 

kept  the  same,  at  53MPa,  during  the  entire  experiment.  The  viscosity  change 

as  a  result  of  structural  relaxation  during  the  course  of  a  viscosity 

measurement  can  be  calculated  from  qigo  and  the  relaxation  rate,  n>  determined 

from  isothermal  annealing  experiments  [10].  For  example,  during  the  30 

minutes  test  time  to  measure  point  27,  the  viscosity  is  calculated  to  have 

changed  by  0.2%.  Optical  microscopy  of  the  tested  specimen  revealed  <10  vol .% 

crystallinity,  originating  from  the  specimen  surface.  This  amount  of 

crystallinity  is  known  to  increase  the  viscosity  gradually  by  at  most  ~25% 

[10]  and  to  have  very  little  effect  on  the  magnitude  of  Q. 

iso 

4.  DISCUSSION 

Three  basic  observations  can  be  made  about  the  isoconf igurational  data 
shown  in  Figs.  1-3:  (i)  the  temperature  dependence  of  isoconf igurational  flow 
is  inherently  non-Arrhenian,  (ii)  Q,gQ(T)  is  dependent  on  the  thermal  history, 
and  (iii)  the  isoconf igurational  viscosity  curves  exhibit  so-called  crossover 
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behavior.  The  implications  of  these  observations  are  discussed  below. 

The  experimental  isoconf igurational  data  for  q(T)  are  non-Arrhenian  and 

therefore  cannot  be  interpreted  simply  in  terms  of  the  temperature  dependence 

of  the  jump  frequency  k^1  (see  equation  (1)).  The  temperature  dependence  of 

the  equilibrium  viscosity,  n  ( T) .  can  be  described  well  by  the  temperature 

eq 

dependence  of  n,(T)  only  [10,11];  thus  the  temperature  dependence  of  k^'  is 
probably  quite  weak.  The  isoconf igurational  activation  energy  must  therefore 
be  accounted  for  by  changes  in  n^  and  hence  by  some  structural  changes  between 
each  of  the  tests  of  an  isoconf igurational  set.  Since  the  viscosity  values 
can  be  reproduced  upon  temperature  cycling,  these  structural  changes  must  also 
be  reversible. 

The  precise  microscopic  picture  which  leads  to  reversible  changes  in  n^ 
is  still  unclear.  It  appears,  however,  that  some  redistribution  of  the  free 
volume  change,  resulting  from  the  thermal  expansion  accompanying  a  test 
temperature  change,  must  occur  [7].  Microscopically,  the  structural  changes 
could  be  described  as  fast,  slight  adjustments  in  the  atom  positions,  without 
changes  in  nearest  neighbors  or  chemical  order. 

The  isoconf igurational  data,  particularly  those  shown  in  Fig.  3,  clearly 

illustrate  the  thermal  history  dependence  of  Qis0(T).  At  testing  temperatures 

far  from  T^  (T  £  560K),  Qjl0  was  lower  if  the  fictive  temperature  of  the 

specimen  was  lower,  i.e.  if  the  specimen  was  more  ordered.  For  example, 

around  5351  the  magnitude  of  Q^#0  ranged  between  0.4  and  2.0eV,  depending  on 

the  state  of  annealing.  For  example,  after  annealing  the  same  specimen  above 

its  fictive  temperature  (after  set  C  in  Fig.  3),  thereby  creating  more 

disorder,  the  isoconf igurational  curve  D  was  obtained  which  exhibited  a  higher 

Qi#0(T).  Above  5  80K,  0 remained  nearly  constant  at  3.6  +  .3eV,  regardless 

of  thermal  treatment.  A  temperature  dependent  Q.  has  also  been  found  in 

iso 
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amorphous  Mg2Si18  [7]  . 

The  history  independent  Arrhenian  behavior  observed  by  other 
investigators  [3,4,5 ,6]  is  reconciled  with  the  results  reported  here  as 
follows.  The  fictive  temperatures  of  the  specimens  in  the  earlier  tests  were 
likely  still  quite  high,  since  the  preannealing  treatments  were  performed  at 
temperatures  well  below  T  .  For  example,  Taub  and  Spaepen  [4]  preannealed 
their  Pd,^  jCu^Si^  j  specimen  at  temperatures  between  500  and  552K  for 
varying  lengths  of  time,  totalling  about  24  days.  Therefore,  their  data  are 
expected  to  exhibit  only  a  slight  deviation  from  apparent  Arrhenian  behavior 
which  is  not  likely  to  be  noticed  over  their  temperature  span  of  40K  or  less. 
Note  that  their  data  (shown  as  the  dashed  curve  from  520  to  550K  in  Fig.  3) 
have  a  temperature  dependence  similar  to  those  shown  as  curve  B  in  Fig.  3, 
which  were  obtained  from  a  specimen  preannealed  for  a  much  shorter  time,  but 
at  higher  temperatures. 

In  the  free  volume  model  [12],  neq(T)  is  determined  by  the  average  free 

volume,  v^.  Crossover  of  isoconf igurational  curves  obtained  from  the  same 

specimen  after  various  degrees  of  annealing  treatments,  for  example  as  in  the 

data  shown  in  Figs.  1  or  3 ,  implies  that  the  same  number  of  flow  defects  can 

be  present  in  specimens  of  different  average  free  volnes.  Therefore,  q.  (T) 

iso 

cannot  be  a  function  of  only  a  single  structural  parameter,  such  as  v^,  but 
that  a  second  parameter  such  as  the  shape  of  the  free  volume  (nonequilibrium) 
distribution  curve  must  be  taken  into  account. 


5 .  CONCLUSIONS 

The  experimental  results  summarized  in  this  paper  show  that 
isoconf igurational  flow  is  inherently  non- Arrhenian.  Hence,  the  apparent 


wW*r a.*  V*  O  •_*  X-*. 
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activation  energy  Qi$0(T)  cannot  *  imply  be  ascribed  to  the  temperature 
dependence  of  the  jnap  frequency  kf*  and  must  be  accounted  for  by  changes  in 
between  each  of  the  tests  in  an  isoconf igurational  set.  Since  ti£#0(T)  i* 
reproducible  upon  temperature  cycling,  the  changes  in  n^  must  also  be 
reversible.  **  thermal  history  dependent.  The  value  of  Q£#0(T)  at 

T  2  5 80S  was  found  to  remain  nearly  the  same  (3.6  +  .3eV)  regardless  of 
annealing  treatment,  whereas  Qijo(T)  at  T  i  560K  decreased  with  increasing 
degree  of  order  in  the  specimen. 

Crossover  of  isoconf igurational  curves  obtained  from  the  same  specimen  at 
various  stages  of  annealing  imply  that  nij0(T)  cannot  be  a  function  of  only  a 
single  structural  parameter. 
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FIGURE  CAPTIONS 


Figure  1.  Isoconf igurational  viscosity  measurements  on  s  Pd__  cCu*Si,*  . 

ribbon  preannealed  under  load  at  S29E  for  457  hours  prior  to  the 
measurement  at  point  1.  The  numbers  indicate  the  order  of  the 
tests.  The  light  solid  line  is  the  equilibrium  viscosity  curve 
[10]  . 

Figure  2.  Isoconf igurational  viscosity  measurements  on  a  Pd ^  gCugSijg  5  wire 
preannealed  under  load  at  593E  for  73  hours  prior  to  the 
measurement  at  point  1. 

Figure  3.  Isoconf igurational  viscosity  measurements  on  a  Pd??  jCu^Sijg  § 
ribbon  preannealed  under  load  at  5931  for  4  hours  prior  to  the 
measurement  at  point  1.  The  dashed  line  represents  earlier 
measurements  by  Taub  and  Spaepen  [4].  The  light  solid  line  is  the 
equilibrium  viscosity  curve  [10]. 
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Figure  1 

Isoconf igurational  viscosity  measurements  on  a  Pd77  5Cu6Silfi  «.  ribbon  preannealed 
under  load  at  529K  for  457  hours  prior  to  the  measurement  at  point  1.  The  nuabers 
indicate  the  order  of  the  tests.  The  light  solid  line  is  the  equilibrium 
viscosity  curve  (10). 
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Figure  3 


Isoconf igurational  viscosity  neasurenents  on  a  Pd77  jCu6Si16  .  ribbon  preannealed 
under  load  at  593K  for  4  hours  prior  to  the  measurement  at  point  1.  The  dashed 
line  represents  earlier  measurements  by  Taub  and  Spaepen  [4].  The  light  solid  line 
is  the  equilibrium  viscosity  curve  [10], 
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ABSTRACT 


The  non-linear  viscosity  change  as  a  result  of  structural 
relaxation  near  equilibrium,  at  temperatures  below  and  above  the 
fictive  temperature,  has  been  measured  for  amorphous 
Pd n  jCugSijg  The  viscosity  increase  is  best  described  by 

bimolecular  annihilation  kinetics  of  flow  defects.  The  viscosity 
decrease  is  described  equally  well  by  either  unimolecular  or 
bimolecular  creation  kinetics  of  flow  defects. 

From  the  asymptotic  values  of  the  model  fits  to  the  data,  the 
equilibrium  viscosity  was  determined.  Over  the  temperature  range 
608  to  630E,  the  temperature  dependence  of  the  equilibrium 
viscosity,  q  ,  of  Pd,^  jCu^Si^  5  **  described  by  the 

Fulcher- Vogel  equation. 


"Present  address:  Sandia  National  Laboratories,  Albuquerque,  NM  87115,  U.S.A. 
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1 .  INTRODUCTION 


Ao  as-prepared  metallic  glass  is  thermodynamically  unstable:  it  can 
lover  its  free  energy  continuously  by  a  series  of  structural  changes  toward 
the  metastable  equilibrium  configuration.  This  phenomenon  is  called 
''structural  relaxation* ' ,  and  is  manifested  by  continuous  changes  in  all 
physical  properties.  The  largest  changes  are  observed  for  the  atomic 
transport  coefficients  such  as  the  shear  viscosity.  q,  and  atomic  diffuaivity, 
D.  Most  measurements  of  the  effect  of  relaxation  on  the  viscosity  have  so  far 
been  made  at  temperatures  veil  below  the  glass  transition  temperature  to  avoid 
crystallization.  At  these  low  temperatures,  where  the  glass  is  necessarily 
far  from  equilibrium,  the  viscosity  has  been  observed  to  increase  linearly 
with  time,  i.e.  q  is  a  constant  [1,2,3,41.  This  is  illustrated  schematically 
by  line  AB  in  Fig.  1.  However,  as  demonstrated  by  Chen  and  Turnbull  for  a 
AuGeSi  metallic  glass  [5],  when  the  structure  of  the  glass  approaches  its 
equilibrium  configuration  (line  DE  on  Fig.  1),  the  rate  of  the  viscosity  rise 
decreases,  since  the  viscosity  must  arrive  at  a  unique,  finite,  equilibrium 
value,  q^(T)  (point  E  on  Fig.  1).  Although  measurements  of  itself  have 
been  reported  for  metallic  glasses  (for  example  on  AuGeSi),  prior  to  our  work 
described,  no  quantitative  experimental  study  has  been  made  of  the  viscosity 
relaxation  kinetics  of  the  approach  to  q#^  near  equilibrium. 

In  this  paper,  we  present  new  measurements  of  viscosity  relaxation  near 

equilibrium,  and  compare  the  results  with  the  predictions  of  a  number  of 

phenomenological  models.  Measurements  near  equilibrium  were  made  possible  by 

using  a  crystallization-resistant  amorphous  alloy,  Pd_,  ,Cu*Si, ,  ..  Non- 

77  *3  O  10  (S 

e 

linear  viscosity  relaxation  (i.e.,  q  not  constant  in  time)  from  both  below  and 
above  q#^  was  observed  (respectively  DE  and  CE  in  Fig.  1).  The 
phenomenological  description  of  the  relaxation  near  equilibrioa  is  shown  to  be 
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in  accord  with  the  linear  viscosity  increase  observed  far  froai  equilibrium. 

The  functional  form  of  n  (T)  is  of  fundamental  interest  as  a  test  of 

eq 

theoretical  models  of  liquids  and  glasses.  Because  many  metallic  glasses  tend 

to  crystallize  before  reaching  their  fully  relaxed  state,  the  equilibrium 

viscosity,  n  (T),  is  often  very  difficult  to  measure  directly.  Ve  will 
eq 

demonstrate  how  values  of  q  (T)  can  be  extracted  from  the  extrapolated 

eq 

asymptotic  behavior  of  the  relaxation  data.  Coincidence  of  the  asymptotic 
values  for  increasing  and  decreasing  viscosities  (DE  and  CE  on  Fig.  1) 
provides  an  experimental  check  on  the  accuracy  of  the  determination  of  n 

eq 

A  preliminary  report  on  the  early  stages  of  this  work  has  been  published 
[6].  In  the  present  paper,  the  complete  data  and  the  full  quantitative 
analysis  of  both  the  relaxation  kinetics  and  equilibrium  viscosity  are 
presented. 


2 .  EXPERIMENTS 

The  material  used  in  these  experiments  was  glassy  Pd ^  jCugSi^g  g 
produced  by  disk  quenching  as  a  long,  continuous  ribbon,  13|im  thick  and  0.22mm 
wide,  from  which  all  specimens  were  cut.  The  amorphous  nature  of  the 
specimens  was  verified  by  x-ray  diffraction  and  differential  scanning 
calorimetry.  The  glass  transition  and  kinetic  crystallization  temperature 
were  measured  to  be  636  and  674E,  respectively  at  20K/min. 

The  viscosity  measurements  were  made  using  the  tensile  creep  apparatus 
described  previously  [7].  Due  to  the  high  test  temperatures,  the  tests  were 
performed  in  air;  no  silicone  oil  was  used.  Nevertheless,  the  spatial 
temperature  profile  of  the  furnace  remained  very  satisfactory:  the 
temperature  variation  within  the  defined  gage  length  of  3.1cm  was  only  £.1/41, 
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and,  although  the  spatial  teaiperature  profile  was  not  an  ideal  step  function, 
the  strain  contribution  from  the  specimen  length  outside  the  gage  length  was 
negligible,  (less  than  1%  of  the  total  strain),  due  to  the  strong  temperature 
dependence  of  the  strain  rate.  The  warm-up  time,  from  the  holding  temperature 
of  314K  to  the  testing  temperature,  varied  from  5  to  13  minutes,  depending  on 
the  testing  temperature  and/or  the  controller  settings.  The  temperature  was 
maintained  to  +1/4K  during  the  anneals,  which  lasted  from  a  few  minutes  to 
about  a  day.  The  furnace  temperature  was  calibrated  with  the  melting  point  of 
pure  lead. 

The  high  temperature  runs  typically  produced  total  creep  strains  of  3-5%. 
The  corresponding  stress  increase  due  to  the  cross-sectional  area  decrease  was 
partially  offset  by  geometrical  changes  in  the  creep  apparatus.  The  error 
introduced  by  the  assumption  of  constant  stress  is  therefore  quite  small.  The 
equivalent  shear  stress  ranged  from  7  to  35  MPa.  Our  experience  showed  that 
this  stress  range  was  still  well  within  the  Newtonian  viscous  regime. 

Loading/unloading  experiments  were  performed  to  ascertain  that  the 
observed  viscosity  behavior  during  the  isothermal  anneal  is  not  affected  by 
changes  in  the  stress-strain  rate  sensitivity,  dlny/dlnr,  the  stress 
dependence  of  q  or  anelastic  transient  decay. 

The  effects  of  stress  changes  on  the  viscosity  of  a  specimen  at  621K  is 
shown  in  Fig.  2.  Note  that:  (1)  at  this  temperature  the  anelastic  strain 
completely  decays  in  a  few  minutes,  (2)  at  these  levels  of  stress  and  stress 
changes,  the  total  anelastic  strain  is  very  small  compared  to  the  plastic 
strain,  and  (3)  the  relaxation  kinetics  are  stress  independent  over  the  range 
of  the  test;  save  for  short-lived  deviations  due  to  the  presence  of  the 
transient  after  a  load  change,  a  smooth  curve  can  be  traced  along  the  curves 
for  q  at  different  stress  levels. 
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The  incremental  anelastic  strain  Ay^  is  determined  by  extrapolation  of 

the  steady-state  strain  bach  to  the  time  when  the  stress  change.  At  was  made. 

—11  2 

An  anelastic  compliance.  Ay^/At  =  8.7  +  2.7  x  10  m  /N  is  obtained  from  the 
results  in  Fig.  2,  in  good  agreement  with  Chen  and  Goldstein’s  measurements 
for  the  same  alloy  [8].  Results  of  tests  at  other  temperatures  are  summarized 
in  Table  1 . 

The  ' 'noise'*  on  the  viscosity  curve  of  Fig.  2  is  a  result  of  the 
differentiation  procedure  used  to  determine  the  slope  of  the  elongation  curve. 
The  continuous  decrease  of  this  slope  leads  to  the  increase  in  the  signal-to- 
noise  ratio  with  time,  seen  on  Fig.  2  which  is  an  experimental  artifact 
without  physical  significance. 


3.  KINETICS  OF  VISCOSITY  RELAXATION 


3-1  Defect  Relaxation  Models 

In  the  Newtonian  regime,  the  viscosity  can  be  written  as  [9] : 


(1) 


( nf :  concentration  of  microscopic  sites  producing  shear,  or  ’’flow  defects’’; 

kf';  stress  unbiased  jump  frequency  of  a  shear  defect;  y0*.  shear  strain  per 

jump  of  a  defect:  vo:  volume  of  a  defect).  Theories,  such  as  the  free  volume 

or  configurational  entropy  models,  can  account  for  q  (T)  over  a  very  large 

eq 

viscosity  range  entirely  by  a  change  in  the  equilibrium  flow  defect 
concentration  n^(T) .  It  is  therefore  reasonable  to  describe  the 
nonequilibrium  viscosity  changes  associated  with  structural  relaxation  by  the 


annihilation  (or  creation)  of  excess  (or  deficit)  flow  defects.  Annihilation 

of  flow  defects  can  be  thought  to  occnr  at  special  sites  called  relaxation 

defects.  The  rate  of  disappearance  of  flow  defects,  n^,  is  then  proportional 

to  the  product  of  the  concentration  of  flow  and  relaxation  defects,  n^nr. 

However,  since  the  system  in  equilibrium  has  a  finite  viscosity,  q  ,  a  number 

eq 

of  defects,  n#  ,  must  remain  after  an  infinite  annealing  time.  Kinetic  laws 

X  s 

in  which  n^  becomes  zero  when  equilibrium  is  reached  (n£  =  )  can  be 

formulated,  such  as  first-  and  second-order  rate  equations  in  the  excess 
number  of  defects  (n£  -  »£>e_).  Respectively,  for  the  unimolecular  model: 


d(nf~nf.ec 

dt 


‘  kl(nrnf.eq> 


and  for  the  bimolecular  model: 


d(nf-nf  ) 

_ £ _ Lt«a_ 

dt 


"  k2(“r»f,eq>‘ 


For  the  decreasing  viscosity  data,  k2  is  negative.  Equation  (2)  describes  an 

annihilation  mechanism  where  n^  is  constant  whereas  equation  (3)  describes  an 

annihilation  mechanism  where  n  <*  (n,-n,  ).  If  n  *  n, ,  a  third 

r  i  i , eq  r  t 

annihilation  equation  is: 


d(nf-nf  .„) 
_ 1 _ f  >  cq 


-  k3nf(nf-n£>eq) 


Integrating  equations  (2),  (3)  and  (4)  gives  for  the  viscosity,  respectively: 
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(6) 


n  =  ^eq  *  <T»eq  "  n0)exp(-k3t)  (7) 

Note  that  when  q  >>  q(t),  as  in  low  temperature  relaxation  experiments,  or 
when  k^  is  small  enough  to  allow  first-order  expansion  of  the  exponential, 
all  the  above  equations  predict  a  linear  viscosity  increase  with  time. 

The  data  exhibiting  non-linear  viscosity  relaxation  were  fit  mostly  to 
the  unimolecular  and  bimolecular  equations  (eqs.  (5)  and  (6))  for  two  reasons. 
First,  the  bimolecular  and  unimolecular  models  have  been  shown  to  describe  the 
viscosity  increase  and  decrease,  respectively,  in  silicate  glasses  [10,11]. 

The  physical  basis  for  this  is  discussed  further  in  Section  S.  Second,  for 
several  sets  of  increasing  viscosity  data  fit  to  all  three  annihilation 
equations,  the  quality  of  the  fit  to  eq.  (7)  obtained  was  intermediate  between 
that  of  the  other  two  equations.  Therefore,  fitting  all  the  data  sets  to  the 
third  equation  judged  unlikely  to  give  significant  new  information. 


Since  a  quantitative  comparison  of  the  unimolecular  and  bimolecular  fits 
is  difficult  with  the  noisy  differentiated  experimental  viscosity  curve,  the 
fits  were  performed  on  the  original  elongation  data.  Equations  (5)  and  (6) 
were  therefore  integrated  to  give  expressions  for  the  elongation  change  with 
time,  respectively: 


Xo  +  * 


%qVl 


(l-exp( 


-  v»] 
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l(t)  =  10  +  I 


[(T»eq'l»o) 

"2.,‘ 


(9) 


K  is  the  product  of  the  sheer  stress  end  specimen  gege  length;  b  is  equel  to 
(ne<1-n0)k2^,,eq*  ^ree  *ree  parameters,  q  ,  or  ^  *n^  the  initiel 
viscosity,  no.  together  with  the  fixed  experimentel  velues  for  1Q  end  K  were 
used  in  fitting  the  dete.  A  non- linear,  multi-parameter  fitting  algorithm 
developed  by  Marquardt  [12]  was  used.  Although  qQ  can  be  obtained  from  the 
data,  it  was  treated  as  a  fitting  parameter  to  avoid  biasing  the  fits  with  any 
transients  which  could  have  been  present  at  the  beginning  of  the  anneal.  If 
the  initial  viscosity  obtained  from  the  fit  was  too  different  from  the 
measured  one,  indicating  the  presence  of  transients,  the  first  few  points  of 
the  run  were  deleted.  Refitting  the  shortened  data  set  then  results  in  closer 
agreement  between  the  fit  and  experimental  values  for  t|o> 

Figure  3  shows  the  viscosity  increase  in  a  Pd ^  jCugSi16  s  as-quenched 
ribbon,  during  isothermal  annealing  at  529K.  For  the  entire  length  of  the 
test,  the  plot  is  linear;  it  shows  no  evidence  of  saturation  to  a  constant 

e 

value.  The  temperature  dependence  of  the  slopes  of  these  plots,  q(T),  is 
shown  on  the  low  temperature  side  of  Fig.  4. 

Figure  5  shows  the  viscosity  increase  for  an  as-quenched  PdCuSi  ribbon 
specimen  at  615E.  The  viscosity  increase  is  clearly  non-linear,  with  q  first 
decreasing  as  the  equilibrium  state  is  approached,  and  then  increasing  due  to 
the  onset  of  crystallization. 

To  determine  the  point  at  which  crystallization  effects  were  no  longer 
negligible,  an  increasing  fraction  of  each  annealing  data  set  was  fit  to  the 
relaxation  models.  In  a  perfect  isothermal  annealing  experiment,  i.e.,  one  in 
which  the  sample  remained  fully  amorphous,  the  same  set  of  fitting  parameters 


/V.V.V.V.V.,..- 


-25- 


should  be  obtained  regardless  of  the  fraction  of  the  test  used.  The  point  at 

which  crystallization  effects  become  significant  is  taken  to  be  that  at  which 

u  end  begin  to  increase  sharply  from  a  reasonably  constant  walne  (see 
eq 

Fig.  6) .  X^  is  defined  as  the  sum  of  the  squared  difference  between  the 

calculated  and  experimental  elongation  data. 

In  all  the  increasing  viscosity  experiments  near  Tg,  the  X*  of  the 

unimolecular  fit  was  larger  than  that  of  the  bimolecular  one.  Unlike  the 

bimolecnlar  *  which  exhibited  a  reasonably  constant  value  before  rapidly 

increasing  due  to  incipient  crystallization,  the  unimolecular  x*  increased 

with  increasing  fraction  of  the  data  included  in  the  fit.  Figure  6  shows  a 
2 

plot  of  t)  and  X  lor  the  unimolecular  and  bimolecular  fits  as  a  function  of 
eq 

the  longest  annealing  time  included  in  the  fit,  for  the  specimen  in  Fig.  5. 

For  this  specimen,  appreciable  crystallization  was  judged  to  occur  after  about 
4  hours,  based  on  the  rapid  departure  from  the  reasonably  constant  initial 
value  of  neq.  The  value  of  i\e^  plotted  in  Fig.  14  was  averaged  from  the  qeq 
value  between  2  and  4  hours  of  the  run.  It  was  assumed  that  for  a  given  run 
the  scatter  in  n  ,  with  increasing  fraction  of  the  data  (with  no  appreciable 
crystallization)  included  in  the  fit,  is  due  to  random  errors.  The  error  bar 
indicated  in  Fig.  14  was  calculated  as  (ti  “*x  -  n  “in)/4,  and  is  assumed  to 
be  the  standard  deviation.  From  crystallization  studies  (discussed  below),  a 
few  percent  crystallinity  is  known  to  be  present  at  3  hours  into  the  anneal. 
Figure  5  shows  the  fits  obtained  from  the  unimolecular  and  bimolecular  models 
as  light  and  heavy  curves,  respectively.  Only  the  first  three  hours  of  the 
creep  test  were  used  for  fitting.  The  bimolecular  model  fits  the  data  very 
well  whereas  the  unimolecular  model  does  not,  due  to  the  fundamental 
difference  in  shape  between  the  model  viscosity  curves.  Figure  7  is  an 
enlarged  plot  of  the  first  three  hours  of  Fig.  5.  Figure  8  shows  the  model 


-26- 

2 

fits  to  the  original  elongation  data.  For  the  fits  shown,  the  unimolecular 
is  17  tiaes  larger  than  the  biaolecalar  x2 •  Figure  9  shows  the  viscosity 
increase  for  another  specimen  at  620K,  again  illustrating  the  superiority  of 
the  biaolecular  fit. 

We  have  also  observed  the  isothensal  approach  to  equilibrium  in  the 

reverse  direction:  a  speciaen.  annealed  for  a  long  time  at  a  low  temperature 

(A  ->  B  on  the  schematic  of  Fig.  1),  was  brought  rapidly  to  a  temperature 

beyond  the  equilibrium  curve  (B  -)  C  os  Fig.  1)  and  the  viscosity  decrease  was 

observed  (C  *>  E  on  Fig.  1) .  Figure  10  shows  the  viscosity  decreasing  for  a 

specimen  first  annealed  at  S91E  for  2  days  before  being  brought  to  615K  at  a 

heating  rate  of  24K/min.  The  viscosity  eventually  increases  due  to 

crystallization.  The  data  before  crystallization  can  be  fit  very  well  using 

either  the  unimolecular  or  bimolecular  model.  For  the  first  0.64  hours  of  the 

2  2 

data,  the  unimolecular  X  is  1.01  times  the  bimolecular  x  •  Both  fits  give 
nearly  the  same  asymptotic  value  for  qeq  that  is  also  in  excellent  agreement 
with  the  asymptotic  value  obtained  from  the  bimolecular  fit  to  increasing 
viscosity,  as  shown  in  Figs.  11  and  12.  Note  that  the  apparent  viscosity 
decrease  is  faster  for  the  decreasing  viscosity  case.  The  same  difference  has 
been  observed  by  Lillie  in  his  viscosity  aeasureaents  on  soda-lime  glass  [10]. 

Figure  13  shows  the  viscosity  behavior  for  a  speciaen  first  annealed  at 
594K  for  164  hours  before  being  brought  to  608E.  The  apparent  viscosity 
remained  constant  before  eventually  increasing  due  to  crystallization.  The 
constant  viscosity  value  was  taken  to  be  q#^. 

The  biaolecular  rate  constants,  k2>  obtained  froa  fits  to  the  increasing 
viscosity  data  are  plotted  in  Fig.  4  together  with  q  (constant)  values 
obtained  by  creep  and  stress  relaxation  aeasureaents  [3].  For  q#^  >>  q(t), 

*2  “  q •  Note  that  k2  has  a  non-Arrheni an  teaperature  dependence.  This 
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behavior  has  also  been  observed  in  Fe4oNi40P14B6  by  Ttnb  tnd  L«b<>**ky  [3]  . 

3.3  Effect  of  Crystallization 

The  viscosity  change  as  a  result  of  structural  relaxation  and 

crystallization  was  fit  to  the  Einstein  equation  for  the  viscosity,  of  a 

mixture  consisting  of  a  small  volume  fraction,  f  ,  of  spherical  particles 

c 

suepended  in  a  medium  of  viscosity  it  [13]: 

^ef f  *  nU  +  2.5fc)  (10) 

The  crystals  were  assumed  to  grow  spherically  with  a  constant  growth  rate,  v, 
from  randomly  positioned  nuclei.  The  Johnson-Mehl-Avrami  analysis  [14]  for 
isothermal  nucleation  and  growth  then  predicts: 

f  (t)  *  1  -  exp(-  at3)  (ID 

where  a  =  N(4/3)nv  ;  N  =  concentration  of  nuclei. 

The  dashed  curves  in  Figs.  5  and  10  are  the  fits  to  the  data,  treating  a 

as  a  new  fitting  parameter  in  combination  with  the  previously  found 

bimolecular  values  for  rj  ,  n  and  k. .  The  calculated  crystalline  volume 

o  * 

fraction  at  the  end  of  each  of  the  tests  in  Figs.  5  and  10  was  0.46  and  0.13 
respectively,  in  reasonable  agreement  with  direct  observation  by 
metallography. 

Metallography  also  revealed  that  crystallization  in  theae  specimens 
originated  from  the  specimen  surface.  Bowever,  a  model  baaed  on  planar  growth 
of  the  crystalline  phase  and  the  same  strain  rates  for  the  amorphous  and 
crystalline  phases  cannot  adequately  fit  the  data  [15,16].  The  remarkable  fit 
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of  equation  10  despite  the  non- random  arrangement  of  non- spherical  crystals, 

is  probably  dne  to  effectively  treating  the  geometrical  factor  2.5  as  part  of 

a  fit  parameter,  since  it  appears  as  a  product  with  f  . 

c 

Patterson  and  Jones  [17]  also  found  excellent  agreement  between  the  creep 

rate  measured  during  crystallization  of  Fe4oNi40P14B6  *nd  tho,e  predicted  by 

eq.  10  up  to  at  least  30%  crystalline  volume  fraction.  In  their  study,  the 

specimen  was  assumed  to  be  at  equilibrium  with  q  eaual  to  the  initial 

eq 

viscosity.  The  change  in  the  crystalline  volume  fraction  was  calculated  from 
the  calorimetrically  measured  crystallization  exotherm  of  the  same  material. 


4.  TEMPERATURE  DEPENDENCE  OF  THE  METASTABLE  EQUILIBRIUM  VISCOSITY 

Figure  14  summarizes  the  neq  values  from  608K  to  630E  which  were  obtained 

in  either  of  three  ways:  (1}  from  the  bimolecular  fits  to  the  increasing 

viscosity  data,  (2)  from  the  averaged  q  values  from  the  bimolecular  and 

eq 

unimolecular  fits  to  decreasing  viscosity  data,  and  (3)  directly  measured 
values  from  experiments  of  the  type  shown  in  Fig.  12.  The  results  of  the  fits 
to  an  Arrhenius  equation  q  •  A  exp(B/T),  Fulcher-Vogel  equation 
q  =  A  exp(B/ (T-Tq) ) ,  and  ''hybrid*'  equation  q  ■  A  exp(B/ (T-TQ) )  exp(C/T)  are 
also  shown.  The  Fulcher-Vogel  and  ''hybrid''  equations  describe  the  data 
well,  whereas  the  Arrhenius  equation  does  not.  The  best-fit  Fulcher-Vogel 
curve  has  A  ■  (1.20  ±  0.4)xl0®  Nsm~2,  B  *  338  +  20K,  Tq  -  583.4  +  IK.  The 
additional  parameter  introduced  in  the  ''hybrid'*  equation  does  not  improve 
the  fit  significantly:  its  •  normalized  for  the  number  of  fit  parameters,  is 
higher  than  that  of  the  Fulcher-Vogel  equation.  Also  shown  in  Fig.  14  are 
earlier  measurements  of  the  equilibrium  viscosity  reported  by  Chen  and 
Goldstein  [81  for  the  same  alloy. 


5.  DISCUSSION 


The  creep  elongation  data  were  fit  to  predictions  of  uniaolecular  and 
bimolecnlar  annihilation  or  creation  of  defects  responsible  for  flow.  It  was 
fonnd  that  the  bimolecnlar  model  described  the  viscosity  increase  dnring 
relaxation  best. 

One  aay  ask  whether  first-  or  second-order  rate  equations  in  the 
viscosity  departure  from  equilibrium: 

dt  “  "  VmV  or  2  (12) 

Bight  not  describe  the  data  just  as  well.  It  turns  out  that  a  second-order 
equation  in  (q^—n)  results  in  an  equation  identical  to  the  biaolecular 
defect-based  one,  whereas  a  first-order  equation  in  results  in  an 

equation  identical  to  equation  (7). 

In  amorphous  covalently  bound  systems,  such  as  soda-lime  glass  [10],  the 
viscosity  increase  during  relaxation  has  also  been  shown  to  be  described  best 
by  the  bimolecnlar  model.  Density  and  refractive  index  studies  above  the 
fictive  temperature  have  shown  essentially  unimolecular  kinetics,  with 
deviations  in  the  early  states  [11] .  In  these  systems,  which  are  modeled  as 
continuous  random  networks,  the  flow  defects  are  easily  identified  as  broken 
bonds  which  allow  shear  rearrangements.  Since  annihilation  of  these  bonds 
occurs  when  two  of  them  come  together  to  form  a  covalent  bond,  the  bimolecnlar 
kinetic  law  of  equation  (3)  is  expected  to  apply  to  the  viscosity  increase. 
Since  foraation  of  a  pair  of  broken  bonds  only  requires  a  single  bond  to  cone 
apart,  the  unimolecular  kinetic  law  of  equation  (2)  is  expected  to  apply  to 
relaxation  above  the  fictive  temperature. 


For  the  metallic  glasses,  one  can  only  speculate  about  the  nature  of  the 
microscopic  mechanism  that  accounts  for  the  bimolecular  kinetics  during  the 
viscosity  increase  [18] .  Since  the  chemical  short  range  order  in  metal- 
metalloid  glasses  such  as  PdCuSi  is  very  strong,  one  can  construct  a  direct 
equivalent  to  the  broken  bond  mechanism  in  covalent  systems.  Diffraction 
evidence  shows  that  the  metalloid  (X)  atoms  are,  ideally,  surrounded  by 
metal(N)  atoms  only  [19].  Direct  analogues  to  broken  covalent  bonds  would  be 
chemical  defects,  i.e.:  X-X  neighbors  and  complementary  excess  H-M  neighbors, 
provided  these  can  be  shown  to  govern  flow. 

It  is  also  possible  to  construct  purely  topological  bimolecular 
annihilation  mechanisms.  In  that  case  a  flow  defect  is  a  cluster  of  atoms 
which  undergoes  a  shear  rearrangement  upon  a  density  or  free  volume 
fluctuation,  resulting  in  a  local  shear  strain  that  is  transferred  elastically 
to  the  specimen  boundary.  The  free  volume  fluctuation,  and  hence  a  flow 
defect,  can  only  be  eliminated  at  a  special  site  or  ''relaxation  defect'', 
which  is  a  cluster  of  atoms,  wherein  as  a  result  of  a  free  volume  fluctuation, 
atoms  rearrange  and  transfer  at  least  part  of  this  free  volume  out  to  the 
specimen  boundary.  Any  rearrangement  with  a  dilatational  component  can  be 
such  a  relaxation  defect  [9].  The  ''p/n  defects''  proposed  by  Egami  et  al. 
[20]  may  be  an  example  of  these  defects,  provided  they  can  be  shown  to  produce 
flow. 

Another  pair  of  topological  defects  can  be  defined  using  Nelson's  model 
for  the  glassy  state,  in  which  a  perfect  tetrahedral  packing  on  the  surface  of 
a  four-dimensional  sphere  is  mapped  into  flat  three-dimensional  space  by  the 
introduction  of  an  array  of  disci inations  [21].  Most  of  these  disclination 
lines  are  long  and  intertwined  and  hence  kinetically  frozen.  Vacancies  or 
interstitials  in  the  four-dimensional  array,  however,  correspond  to  localized 
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disclination  ''babbles''  of  opposite  sign  that  cm  mutually  snaibilste  [22]. 
Again,  if  these  ''babbles''  can  be  shown  to  prodace  local  shear,  their 
annihilation  may  explain  the  bimolecalar  relaxation  kinetics. 

Since  fitting  oar  decreasing  viscosity  data  to  the  bimolecalar  and 
unimolecnlar  equations  coaid  not  distingaish  between  the  two  models, 
speculation  on  the  mechanisms  for  generation  of  flow  defects  would  be 
premature.  Further  measurements  on  an  even  more  stable  glass,  such  as 
Pd40Ni40P2Q.  are  underway  to  resolve  this. 

The  newly  determined  valnes  of  q  are  well  described  by  the  Fulcher- 

Vogel  and  ’ 'hybrid* ’  equations,  bat  not  by  the  Axrhennius  equation.  This 

result  is  expected  in  metallic  glasses  where  the  flow  process  requires  a 

cooperative  motion  involving  a  number  of  atoms.  The  Fulcher-Vogel  equation  is 

also  the  one  predicted  by  the  free  volume  [23]  and  configurational  entropy 

models  [24].  In  contrast,  the  equilibrium  viscosity  behavior  in  systems  such 

as  pure  silica  where  flow  is  largely  due  to  bond-breaking  exhibits  Arrhenian 

behavior.  Van  den  Beukel  and  Radelaar  [25]  advocate  using  the  ''hybrid'’ 

equation  to  describe  the  equilibrium  viscosity  in  metallic  glasses  in  order  to 

account  for  the  rather  high  apparent  activation  energies  observed  in 

isoconf igurational  flow  experiments  [1].  They  pointed  out  that  Chen  and 

Goldstein's  results  [8]  can  be  equally  well  described  by  the  Fulcher-Vogel  and 

' 'hybrid' '  equations.  Moreover,  by  fixing  the  fit  parameter  C  (see  Section  4) 

at  different  values  they  also  showed  that  the  Chen  and  Goldstein  data  can  be 

equally  well  described  by  different  sets  of  ''hybrid'*  fit  parameters.  The 

effect  of  the  additional  fitting  parameter  contained  in  the  ''hybrid'' 

equation  was  to  lower  the  in  the  fit  to  our  data.  However,  the  normalized 
2 

X  for  the  Fulcher-Vogel  fit  remained  the  lowest. 

The  value  of  Tp  for  the  Fulcher-Vogel  fit  to  our  data  is  similar  to  that 


m 
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found  by  Chen  and  Goldstein.  The  paraaeter  B,  however,  is  considerably 
smaller.  (338K  vs.  3180E);  according  to  the  free  volnsie  Model,  with  a  volume 
thermal  expansion  coefficient  of  4.5  x  [26],  our  value  of  B 

corresponds  to  a  critical  free  volume  fluctuation,  v*,  of  only  6%  of  the 
atomic  volume.  The  pre- exponential ,  A,  in  our  fit  is  very  large. 

Extrapolation  of  the  equation  to  higher  temperatures  gives  unphysically  high 
values  for  the  viscosity,  and  is  therefore  unwarranted.  This  is  not  unusual, 
however,  since  it  is  known  that  for  most  amorphous  systems,  metallic  and  non- 
metal  lie  ones,  a  single  set  of  Fnlcher-Vogel  parameters  is  inadequate  to 
describe  the  viscosity  over  the  entire  temperature  range. 

Our  values  of  q  are  clearly  higher  than  those  measured  by  Chen  and 
Goldstein  [8],  by  almost  an  order  of  magnitude.  Since  the  previous  authors 
assumed  that  the  equilibrium  structure  was  attained  shortly  after  the  decay  of 
the  anelastic  transients,  and  since  the  earlier  values  were  not  confirmed  by 
decreasing  viscosity  relaxation,  the  discrepancy  is  probably  due  to 
insufficient  relaxation  in  the  earlier  experiments. 

The  data  for  q#^  plotted  in  Fig.  14  only  span  the  temperature  range  from 
608  to  630K.  Increasing  viscosity  experiments  were  also  performed  outside 
this  range.  At  higher  temperatures,  crystallization  occurred  too  rapidly.  At 
slightly  lower  temperatures,  crystallization  was  delayed,  but  still  occurred 
prior  to  the  decay  of  the  anelastic  transient.  Since  anelastic  transients 
decay  more  rapidly  at  higher  temperatures,  several  tests  were  performed 
whereby  the  specimen  is  first  annealed  at  a  slightly  higher  temperature  than 
the  testing  temperature  for  a  short  period  of  time.  These  attempts  were 
unsuccessful,  however,  in  sufficiently  delaying  crystallization,  probably  due 
to  enhanced  devitrification  at  the  pre-annealing  temperature.  It  should  be 
noted  that  the  q#<^  values  at  the  two  lowest  temperatures  were  obtained  from 


experiment*  that  consisted  of  an  initial  long  pre-anneal  at  a  low  temperature 
prior  to  rapid  warning  to  final  testing  tenperatnre.  This  latter  procedure 
was  effective  probably  because  during  the  pre-anneal  anelastic  transient  decay 
was  completed  with  little  accompanying  crystallization  due  to  the  low 
temperature,  and  high  initial  viscosity  at  the  final  temperature  helped  retard 
the  growth  of  the  crystalline  phase.  At  still  lower  temperatures,  the 
viscosity  increase  was  linear  over  the  entire  duration  of  the  test. 

6 .  CONCLUSIONS 

The  viscosity  increase  as  a  result  of  structural  relaxation  was  observed 

to  be  linear  in  time  at  temperatures  well  below  T  ,  but  non-linear  at 

8 

temperatures  near  T^,  when  the  glass  is  closer  to  equilibrium.  The  viscosity 
increase  could  be  described  best  by  a  mechanism  involving  bimolecular 
annihilation  of  chemical  or  topological  flow  defects.  The  viscosity  decrease 
as  a  result  of  structural  relaxation  from  above  the  equilibrium  value  was  also 
measured  and  found  to  be  fit  adequately  by  both  unimolecular  and  bimolecular 
models  for  creation  of  flow  defects. 

New  values  of  the  equilibrium  viscosity  of  Pd??  5CufiSi16  5  were  obtained 

in  three  ways:  (1)  from  the  asymptotic  value  of  the  bimolecular  viscosity  fit 

to  increasing  viscosity  data,  (2)  from  the  averaged  asymptotic  values  of 

bimolecular  and  unimolecular  viscosity  fits  to  the  viscosity  decreasing  data, 

and  (3)  from  direct  measurements  on  specimens  previously  well-annealed  at  a 

lower  temperature.  Within  experimental  accuracy,  the  same  asymptotic  value 

for  is  obtained  from  fits  to  increasing  and  decreasing  viscosity  data  at 

the  same  temperature.  The  newly  determined  values  of  n  are  about  an  order 

eq 

of  magnitude  higher  than  the  literature  data  for  the  same  alloy.  They  are 


well  described  by  the  Folcher-Vogel  equation. 


The  effect  of  partial  crystallization  on  the  viscosity  coaid  be  described 
by  a  simple  modification  of  the  Einstein  theory  for  flow  of  dispersions. 


ACKNOWLEDGMENTS 

We  thank  K.F.  Kelton  for  the  use  of  his  fitting  program.  This  work  has 
been  supported  by  the  Office  of  Naval  Research,  Contract  N00014-77-C-0002. 


v?/; 


-35- 


FIG  ORE  CAPTIONS 

Figure  1.  Shear  viscosity  q  in  the  various  stability  regimes:  stable 

equilibrium  above  T^,  metastable  equilibrium,  isoconf igurational 
(BC),  structural  relaxation  (AB, DE, CE) . 

Figure  2.  Total  strain  and  apparent  viscosity  changes  as  a  result  of  both 
structural  relaxation  and  delayed  elasticity  after  loading  and 
unloading  of  the  speciaens. 

Figure  3.  Example  of  linear  viscosity  increase  with  tiae  for  Pd,^  ^Cu^Sij^  j 
at  T  =  529K. 

e 

Figure  4.  Arrhenius  plot  of  the  rate  constants  or  q  (constant)  deterained 

from  stress  relaxation  and  creep  experiments. 

Figure  3.  Example  of  the  nonlinear  viscosity  increase  of  a  PdCuSi  ribbon 

specimen  at  615K.  The  fits  to  the  unimolecular  and  bimolecular 

models  are  shown. 

2 

Figure  6.  q  and  X  obtained  by  fitting  the  data  shown  in  Fig.  5  to  the 
eq 

unimolecular  (open  symbols)  and  bimolecular  (solid  symbols)  models 

as  a  function  of  the  longest  annealing  time  included  in  the  fit. 

Figure  7.  Enlarged  plot  of  the  bimolecular  and  unimolecular  fits  to  the  first 

three  hours  of  the  data  in  Figure  5.  The  bimolecular  model  is  a 

better  fit  to  the  data  than  the  unimolecular  model.  The 

fundamental  difference  in  shape  between  the  model  curves  is  clear. 

Figure  8.  The  bimolecular  and  unimolecular  fits  to  the  first  three  hours  of 

the  original  elongation  data  in  Fig.  5.  The  bimolecular  model  fits 

the  data  better  than  the  unimolecular  model. 

Figure  9.  Example  of  the  non-linear  viscosity  increase  for  a  Pd_,  cCu*Si,,  , 

77  #5  o  lo 

ribbon  specimen  at  620K.  The  light  and  heavy  curves  are  the  fits 
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to  the  oniaolecular  and  bioolecolar  models,  respectively. 

Figure  10.  Example  of  the  viscosity  decrease  with  time  during  isothermal 

annealing  of  PdCuSi  at  615E,  above  the  fictive  temperature.  The 
specimen  vas  preannealed  at  5  91E  for  74  hours.  The  light  and 
heavy  curves  are  the  unimolecular  and  bimolecular  fits  to  the 
data,  respectively.  The  dashed  curve  takes  into  account 
crystallization. 

Figure  11.  Viscosity  relaxation  of  two  PdCuSi  specimens  at  612K.  The  upper 
curve  is  for  a  specimen  preannealed  at  595E  for  49  hours;  the 
lower  curve  is  for  an  as-quenched  specimen.  The  fits  to  the 
decreasing  and  increasing  q  give  the  same  asymptotic  value  for 

\q‘ 

Figure  12.  Viscosity  relaxation  of  two  PdCuSi  specimens  at  615E.  The  upper 

curve  is  for  the  specimen  in  Figure  10;  the  lower  curve  is  for  an 

as-quenched  specimen.  The  fits  to  the  decreasing  and  increasing  q 

give  the  same  asymptotic  value  for  q  . 

eq 

Figure  13.  Viscosity  change  with  annealing  time  at  608E  for  a  PdCuSi  ribbon 
specimen  preannealed  at  394E  for  164  hours.  The  viscosity 
remained  constant  for  several  hours  prior  to  specimen 
crystallization. 

Figure  14.  Summary  of  values  determined  from  model  fits  to  increasing 

viscosity  data  (circles),  decreasing  viscosity  data  (dots),  and 

direct  measurements  on  specimens  preannealed  at  low  temperature 

(circled  x) .  The  temperature  dependence  of  q  is  described  well 

eq 

by  the  Fulcher-Vogel  and  ''hybrid''  equations.  The  Chen  and 
Goldstein  results  are  from  ref.  8. 
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Figure  1 


Shear  viscosity  n  in  the  various  stability  regimes:  stable  equilibrium  above  TM, 
metastable  equilibrium,  lsocanfigurational  (BC) ,  structural  relaxation  (AB.DE.CE) 


Arrhenius  plot  of  the  rate  constants  k2  or  f|  (constant)  determined  from  stress  relaxation  and  creep  experiments 


Figure  5 


Example  of  the  nonlinear  viscosity  increase  of  a  PdCuSi  ribbon  specimen  at  615K. 
The  fits  to  the  unlnolecular  and  blmolecular  models  are  shown. 
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Figure  7 

Enlarged  plot  of  the  biaolecular  and  uniaolecular  fits  to  the  first  three  hours  of 
the  data  in  Figure  5.  The  biaolecular  aodel  is  a  better  fit  to  the  data  than  the 
uniaolecular  aodel.  The  fundaaental  difference  in  shape  between  the  aodel  curves 
is  clear. 
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Figure  8 

The  bimolecular  and  unimolecular  fits  to  the  first  three  hours  of  the  original 
elongation  data  in  Figure  5.  The  bimolecular  model  fits  the  data  better  than  the 
unimolecular  model. 
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Figure  9 

Example  of  the  non-linear  viscosity  increase  for  a  Pd^y^Cu^Siig^  ribbon  specimen 
at  620K.  The  light  and  heavy  curves  are  the  fits  to  the  unimolecular  and  bimolecular 
models,  respectively. 


Example  of  the  viscosity  decrease  with  time  during  isothermal  annealing  of  PdCuSl  at  615K,  above  the  fictive  tempera¬ 
ture.  The  specimen  was  preannealed  at  591K  for  74  hours.  The  light  and  heavy  curves  are  the  unimolecular  and 
bimolecular  fits  to  the  data,  respectively.  The  dashed  curve  takes  into  account  crystallization. 
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Figure  11 

Viscosity  relaxation  of  tvo  PdCuSi  specimens  at  612K.  The  upper  curve  is  for  a 
specimen  preannealed  at  595K  for  49  hours;  the  lower  curve  is  for  an  as-quenched 
specimen.  The  fits  to  the  decreasing  and  increasing  f)  give  the  same  asymptotic 
value  for  n  . 

*9 


st2j: 


0 


2 


4  6 
TIME  (hours) 


8 


IQ 


Figure  12 

■P< claea  in  Figure°10°fthe°l PdCuS1  •P«cimens  at  615k  Th 

t#  th* — ssSSSps^ 

n«q* 


A1IS0DSIA 


",  m 

S 

r.“ 

f  ; 

? 

*% 

< 

| 


-51- 


Table  1 

MEASURED  TIME  FOR  AN ELASTIC  DECAT  AND  AN ELASTIC  COMPLIANCE 


P 


Run 

Test  Temperature 

(R) 

Decay  Time 

(min) 

Aya/At 

(10-11m2/N) 

A 

614 

10  +1 

9.0  +  2.7 

B 

614 

4.4+1 

6.3  +1.2 

C 

614 

8.7  +  2 

10.0  +3.6 

D 

615 

6.1+1 

7.9  +1.5 

E 

621 

2.9  +  1 

8.7  +  2.7 

r« 
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